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Abstract. We discuss the concept of an all-optical and ionizing matter- wave 
interferometer in the time domain. The proposed setup aims at testing the wave nature 
of highly massive clusters and molecules, and it will enable new precision experiments 
with a broad class of atoms, using the same laser system. The propagating particles 
are illuminated by three pulses of a standing ultraviolet laser beam, which detaches 
an electron via efficient single photon-absorption. Optical gratings may have periods 
as small as 80 nm, leading to wide diffraction angles for cold atoms and to compact 
setups even for very massive clusters. Accounting for the coherent and the incoherent 
parts of the particle-light interaction, we show that the combined effect of phase and 
amplitude modulation of the matter waves gives rise to a Talbot-Lau-like interference 
effect with a characteristic dependence on the pulse delay time. 
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1. Introduction 

1.1. Interferometry with molecules 

Over the last decades, matter-wave interferometry has developed into a highly active 
research field, ranging from the foundations of physics to quantum enhanced precision 
measurements. Beam splitters, mirrors, diffraction gratings, traps or wave guides are 
nowadays readily available for electrons pQ, neutrons P] and atoms [3], but they still 
represent a substantial challenge for large clusters or molecules. 

Diffraction and interference experiments with molecules were started systematically 
only in the mid-1990s and led to the exploration of far-field diffraction effects 
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at nanomechanical |1H6] and optical gratings [7], as well as at opaque discs [8]. 
Interferometers were operated using nanogratings in a Mach-Zehnder configuration [H], 
running laser waves in a Ramsey-Borde arrangement [IHIITT], mechanical gratings in a 
near-field Talbot-Lau design [T^|T3] or, more recently, in a combination of mechanical 
and optical phase gratings [T^ITS]. 

All these designs have their merits and drawbacks. While far-field experiments are 
conceptually appealing and allow one to spatially separate the diffraction orders, it was 
already pointed out by Clauser [16] that the use of a near-field configuration, such as 
the Talbot-Lau interferometer, offers several advantages, in particular much reduced 
collimation and coherence requirements and therefore increased signals. 

1.2. Optical elements for matter wave interferometry 

Nanomechanical gratings, which serve to block a part of the particle beam, are the most 
natural diffraction elements for matter waves. Using modern lithographic techniques, it 
is nowadays possible to nanofabricate e.g. silicon nitride structures with a precision that 
guarantees a predefined slit period to within a few Angstroms, even over millimeter sized 
areas [TT] . Such masks are of great importance for many applications in atom P^|[TU] and 
electron interferometry [20l|21] since they do not rely on any internal particle property. 

However, for highly polarizable and slow particles, the presence of dispersion forces 
near the grating walls becomes increasingly important |ni[I21[221l23]- These van der 
Waals or Casimir-Polder interactions introduce a phase shift with a strong position 
and velocity dependence [2U[2S]. The force may even remove a substantial part of the 
molecules from the beam when they get too close to the surface. 

In contrast to that, optical standing- wave gratings can neither be destroyed nor 
clogged by large particles. Narrow-band lasers allow one to define the grating period 
with high precision and the grating transmission function, defined by the particle-light 
interaction, can be controlled and modulated in situ and on a short time scale via the 
laser intensity [26|l27]. 

1.3. From phase gratings to absorptive gratings of light 

The use of optical phase gratings is by now well established both in atom interferometry 
(e.g. [27I - I3T] ). with electron beams [32], and also for complex molecules f7\ 111]. 
Interferometry with phase gratings requires a prior effort in preparing sufficient 
transverse coherence in the particle beam. This is naturally provided by an atomic 
Bose-Einstein condensate [33] or if the diffracting particles start already in a periodic 
arrangement predefined by a potential [31]. A spatially extended and incoherent 
molecular beam requires, however, to use an absorptive grating prior to the diffracting 
phase mask to prepare some transverse coherence. An elegant way of realizing an 
absorptive optical grating for atoms is to pump them into a dark state, i.e. into an 
undetected internal state. This strategy works rather well for atoms [55|[5B]. but clusters 
and molecules usually do not offer the needed addressability of selected states. 
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In our present work we focus on the implementation of absorptive gratings via 
single-photon ionization in the anti-nodes of an ultraviolet standing light wave |37j . 
Modern (V)UV lasers readily provide the photon energy required to surpass the 
ionization threshold of many kinds of neutral atoms, clusters and molecules. The optical 
grating will imprint a nanostructured periodic amplitude profile onto the particle beam, 
very much like a material mask and separate charged and neutral particles. We study, 
in particular, the sequence of three ionizing optical grating pulses in the time domain. 
The delay between the pulses serves as a control parameter, while the spatial positions 
of the gratings are fixed. 

The high time resolution of modern lasers allows us to take advantage of all the 
benefits that are shared by other atom interferometers in the time domain [27 t [5n | [55] : 
the interference fringe positions can be largely independent of the particle velocity and 
the interference pattern can be scanned with very high accuracy and without any moving 
mechanical parts. While the idea of ionization gratings is here worked out for massive 
clusters, the method is naturally extended to atoms, as well. 

2. Optical time-domain ionizing matter interferometer 

In the following, we will consider particles that are much smaller than the laser wave 
length, and refer to them as 'clusters', even though the concept also applies to atoms or 
some molecules. The particle-light interaction is then governed by two parameters, the 
(real) optical dipole polarizability a and the absorption cross section o"abs at the laser 
frequency. 

2.1. Interferometric setup 

The layout and time sequence of an Optical Time-domain Ionizing Matter-wave 
(OTIMA) interferometer are depicted in Fig. [1] A pulsed cluster beam passes alongside 
a plane mirror which serves to retro-reflect the beam of an ultraviolet (UV) laser. The 
light source illuminates each cluster cloud in three short pulses, separated by the equal 
time delay T. We denote them by Gi, G2 and G3. The laser beam waist is chosen 
such that all clusters are covered equally by the standing light waves, thus ensuring that 
all particles with the same distance to the mirror experience the same light intensity. 
Depending on the velocity of the cluster beam, this can be realized by a single laser of 
high repetition rate or three pulses that follow the particle cloud along its direction of 
motion. 

The particles are also subjected to a (weak) constant force perpendicular to the 
mirror plane, e.g. gravity or an electric field with VE"^ =const. The field will also 
serve to remove all ions from the beam. Since all relevant forces only depend on the 
cluster-mirror distance, we restrict our attention to this 

A grating pulse acts in several ways: First, the optical dipole potential imprints a 
sinusoidal phase modulation onto the center-of-mass state of the clusters. The period d is 
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Mirror MCP 




Laser gratings 

Figure 1. Experimental layout and pulse sequence for an optical time-domain ionizing 
matter-wave (OTIMA) interferometer. A beam of nanoparticles is chopped into 
bunches, each of which is subsequently illuminated by the same three laser pulses 
(Gi, G2 and G3), separated by a variable but equal delay time T on the order of 
the Talbot time Tx- The retro-reflection of the laser light at a common plane mirror 
generates a phase-stable standing wave. The laser wave length is sufficiently short 
to ionize the particles with high probability at the antinodes after absorbtion of a 
single photon. The ions are removed by a homogeneous electric field. The remaining 
neutral particles fly into a time-of- flight mass spectrometer (TOF-MS), where they are 
photo- ionized, accelerated and detected by a multi-channel plate (MCP). This scheme 
allows one to post-select the clusters with a mass selectivity of better than 0.1 % and 
to perform the experiment with a large number of different masses simultaneously. A 
quantum interference fringe pattern is observed in the count rate by varying the pulse 
separation T, provided the clusters are exposed to an additional constant acceleration 
a, e.g. due to gravitation or a constant electrical field gradient. 



given by half the laser wavelength, c? = Al/2. The maximal phase shift 0o is proportional 
to the real part of the optical polarizability a and to the laser pulse energy. Second, even 
a single absorbed photon will ionize a given cluster. The probability of such ionization 
events is again an oscillating function of the transverse position x. The average number 
uq of absorbed photons in the antinodes is proportional to the absorption cross section 
(Tabs and to the laser pulse energy. The spatially periodic ionization of the clusters, and 
their removal by the electric field, then results in an amplitude modulation of the particle 
beam, much like in a material grating. It is important that a cluster does indeed ionize 
with high probability after the absorption of a single photon. Photoabsorption results 
in a momentum kick of iL/i/Al with a random sign. This could blur the interference 
pattern if the particles were not ionized and removed. A third effect is the radiation 
pressure exerted on the clusters due to Rayleigh scattering. This is a potential source 
of decoherence since the scattered photons impart momentum kicks onto the clusters 
which are comparable to the 'grating momentum' Apc = h/d hut randomly directed (see 
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Sect. I4.ip . Behind the interferometer the remaining neutral clusters are detected mass- 
selectively. The interference pattern will manifest itself in a characteristic dependence of 
the count rate on the particle mass m, the pulse separation T, the external acceleration 
a and also the laser pulse intensity. 



2.2. Near field interference 

Near-field Talbot-Lau interferometers (TLI) are particularly well suited for exploring 
quantum wave mechanics with massive particles and short de Broglie wave lengths 
P ^ IT H ITG I 12 ^ 125] . In our time-domain version, the laser pulses modify the x-component 
of the motional cluster state in the same way as the three gratings of a TLI 'in space'. 
The two concepts are related to each other by a change of the reference frame. 

Compared to a mechanical Talbot-Lau setup, the optical analog adds substantial 
control due to the precisely defined delay between the pulses and the possibility of tuning 
the pulse strengths individually. This allows one to overcome source imperfections, in 
particular the longitudinal velocity spread in the beam, and to optimize the fringe 
visibility. In contrast to the spatial near-field interferometer, which is characterized by 
the Talbot distance Lx = d'^/X, the characteristic scale in the time domain is given by 
the Talbot time [27] 

Tt = md^/h (1) 

This is the typical time for the self-imaging of a pulsed diffraction grating under coherent 
illumination. 

In practice, the x-momentum distribution Ap^ of the particle cloud will be much 
broader than the 'grating momentum' h/d. The role of Gi at t = is therefore reduced 
to acting as a mask that prepares spatial coherence by only transmitting particles from 
many parallel but very narrow 'sources' in the nodes of the standing laser wave. Some 
transverse coherence thus emerges after time t = T, when the second grating pulse G2 
is applied. The amplitude and phase modulation at G2 then leads to a resonant spatial 
modulation of the cluster density at time t = 2T |271|39]. The normalized particle 
density (Sect. [3]) then takes the form 



W2t(X 



This pattern is again periodic in x, with period d but it is shifted in the presence of an 
acceleration a. 

The Fourier components characterizing the fringe pattern ([2]) can be specified in 
terms of the Talbot-Lau coefficients of the k-th grating [15], 



Ccoh Ci( 
Ccoh ~l~ Ci( 



m/2 



Jm (^Sgn(Ccoh + Cion)^Ccoh-Cion^ , (3) 



involving the expressions 

Cion = 4'\os(7rO/2, (4) 



Ccoh = 0o''^ sin «). (5) 
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Here Jm is the Bessel function of the first kind and Uq and 0q are the above mentioned 
grating pulse parameters, specifying the maximal mean number of absorbed photons and 
the maximum phase shift, respectively. They are defined in Sect. Oby Eqs. (ITT!) and 



Note in Eq. ([2]) how the dependence on ^ = £T/Tt in the second grating coefficient 
modulates the parameters (jl]) and ([5]) of ionization and phase shift, and thus relates 
the pulse timing T to the Talbot time Tp. It is in this functional dependence that the 
quantum scale Tp enters, leading to the rich fringe structures characteristic of near field 
interference phenomena. 

The components attributed to the first grating mask [k = 1) take the simpler form 



where Im is the modified Bessel function of the first kind. Note that this factor depends 
only on the photon absorption not on the phase shift. In the limit of small photon 
absorption n^Q^ — )• 0, we find Bm{0) = Sm,o- This implies that if the first optical grating 
does not ionize the clusters the final density distribution will be as uniform as the 
initial particle cloud, even if there is a substantial additional phase shift. In the case of 
equal time delays T between the pulses, the first grating must therefore be ionizing, i.e 
absorptive, in order to generate a fringe pattern. 

Near-field interference may arise behind pure phase gratings, under some 
circumstances: As discussed in Sec. [31 coherent rephasing effects can be expected at 
times T + r after the second pulse, with |r| ^ T. Such transient near- field diffraction 
phenomena were observed in several time-domain experiments [2711101 HI] where the 
Bragg reflection of the third laser pulse was used to detect the reconstructed atomic 
density pattern. In contrast to that, our present proposal uses the symmetric Talbot- 
Lau recurrence with its high visibility and robustness, as required for slow cluster beams 
detected by a third absorptive mask. 

2.3. Fringe visibility 

The fringe pattern is finally probed by the third laser field G3 which ionizes and 
extracts the clusters in its antinodes and transmits the remaining neutrals to the 
detector. Their fraction is proportional to 



When recorded as a function of the transverse acceleration a, this periodic signal has a 
form similar to the density pattern (|2|). The Fourier coefficients now merely contain in 
addition the components Bn\o) of the third grating pulse. A useful interference pattern 
can be recorded even if the force is constant as in the presence of gravity a = g: In this 
case one varies the pulse delay T on the scale AT < d/2gT. If this is much smaller than 




(6) 
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the Talbot time Tt the resulting pattern is practically periodic in since the Fourier 
coefficients then hardly change on that scale. This is naturally the case for massive 
clusters, while one requires a high multiple M ^ 1 of the Talbot time, T = MTt, for 
atoms and light molecules. 

The fringe visibility of the periodic signal ([7]) is conventionally defined as 
q — Q . 

^ = 9 — TV~' ^ ^ 

In practice, the pattern is often close to sinusoidal, allowing one to describe the 
experimental observation by fitting a (i-periodic sine curve with offset So and amplitude 
2Si to the measured data. It is then more appropriate to consider the 'sinusoidal 
visibility' Vsin, the visibility associated to the ffist Fourier component. For the expected 
signal ([7]), and after the insertion of ([3]) and ([6]), it is given by 



h Uo^/2 



(Ccoh ~ Cion)^2 


(VCcoh Cion^ 




ICcoh ~l~ CionI 


lo Uo^/2) 





So I, (4^)/2 ) ICeoh + CionI lo (^'VS ) h ( ^'VS 



(9) 

Note that the sinusoidal visibility may assume values larger than one, while the 
conventional visibility (|8]) by definition cannot exceed 100%. 

2.4- Quantum or classical patterns? 

It is important to keep in mind that the observation of fringes in a two-grating or three- 
grating setup is, by itself, not yet a conclusive proof of quantum interference, since moire- 
type shadow patterns can also be created by particles moving on classical trajectories. 
The genuine quantum origin of the patterns is revealed by their characteristic 
dependence on the quantum scale ([1]) and the detailed functional dependence of the 
fringe visibility on various parameters. It is encoded, formally, in the oscillating in- 
dependence of the Talbot-Lau coefficients ([3]) for the second grating pulse, while the 
actions of Gi and G3 are the same in the quantum and the classical formulation. 
Calculating the deflection of a classical trajectory under the influence of a standing 
light- wave one is led to classical analogs of the Talbot-Lau coefficients [T3]. Their form 
is identical to the quantum coefficients ([3]), except for the definitions @ and ([5]) of the 
parameters Cion and Ccoh- They are now simply given by Cion = ^o^^ /2 and Ccoh = 4>'o^t^^'-, 
in particular, they do not oscillate with ^. 

In agreement with intuitive expectations, the classical and the quantum description 
converge in the limit of large masses, m — )■ 00, or for short pulse separations, T — )■ 0, 
where the argument of the Talbot-Lau coefficients ([3]) tends to zero. They do differ 
markedly for finite masses and times. Genuine quantum interference becomes clearly 
distinguishable from any classical dynamics when the pulse delay time T exceeds at 
least the Talbot time, T > Tt. 
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Figure 2. (a) Simulation of the 'sinusoidal' quantum interference visibility Vsin 

(solid line), the conventional quantum contrast V (dash-dotted line), and the classical 

expectation for the fringe contrast (dashed line) as a function of T/Tt for tiq = 8. 

(2) 

(b) The same contrast curves as a function of Uq for fixed T ^ T^. The simulation 
assumes the material characteristics of gold clusters, /3 = 1.0. 



2.5. Expected visibilities for high-mass clusters 

We now discuss how the expected interference visibihty ([9]) depends on the laser 
pulse parameters and on the optical cluster properties for high-mass metal clusters. 

(k) 

The photon absorption parameters Hq grow linearly with the laser intensities; they 
determine the 'survival factors' bII'\o) for each grating pulse, i.e. the spatially averaged 
probabilities of a particle not to be ionized during the k-th pulse. The optimal laser 
power is therefore a compromise: Large photon numbers may increase the interference 
visibility, but they also decrease the number of clusters that arrive at the detector, 

While the first and the third grating pulse are fully characterized by the absorption 
parameters n^^^ and nl^\ the second pulse shapes the interfering matter wave in 

(2) (2) 

amplitude and phase with parameters Uq and 0o . To account for the specific optical 
properties of the interfering particle independently of the laser intensity it is therefore 
convenient to use the dimensionless ratio 

(2) \ 

P = \ = ^ (10) 

with Al the laser wave length. We assume this to be 157 nm, as will be justified below. 

In order to visualize the expected interference patterns, we keep the first and third 
grating intensities fixed at n^Q ^ = n^^^ = 8 and choose the cluster parameter /3 = 1.0 
to be characteristic for large gold clusters, essentially independently of the particle size. 

Figure [2](a) depicts the sinusoidal visibihty as a function of the delay time T in units of 

(2) 

the Talbot time Tt for Uq = 8. Panel (b), on the other hand, shows the fringe contrast 

(2) (2) 

at T = Tt when the intensity of G2 is varied between Uq = and tZq = 10. Both the 
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Figure 3. Predicted quantum interference contrast as a function of the pulse delay 
time T/Tt for different cluster materials, i.e. for different values of the parameter /3. 
The solid line corresponds to gold (/? = 1.0), the dashed line to cesium (/3 = —1.3) and 
the dashed-dotted line to silver {(3 = 9.2). The /?-values are obtained from Eq. (|17l) , 
using dielectric functions of the bulk materials 42_ . The power of the second pulse is 
set to ni?^ — 8. 



quantum prediction (solid and dash-dotted curves) as well as the classical prediction 
(dashed curves) are included, and the circles in (a) and (b) identify the point of equal 
time and power in both graphs. We observe pronounced visibility peaks recurring at 
multiples of the Talbot time, a feature clearly absent in the classical calculation. In the 
following we therefore consider only the quantum case. Note also that the sinusoidal 
visibility is not always a sufficient approximation to the full fringe contrast (|8]). As seen 
in Fig. 121 it can overestimate the proper fringe amplitude in high-contrast regions (solid 
line), but it can also vanish where the total contrast is still finite (dash-dotted line). 
Depending on the experimental data analysis, higher-order Fourier coefficients may need 
to be accounted for in such cases. 

Figure |3] illustrates the influence of the material parameter /3 on the sinusoidal 
visibility as a function of the pulse separation for three different materials, cesium, gold 
and silver clusters. Their /3- values are representative for most ionizable cluster materials. 
The peaks of the visibility curves resize and move away from integer multiples of the 
Talbot time depending on the magnitude and sign of /3. The general shape, however, is 
mainly determined by the number of absorbed photons Uq and stays roughly the same 
for all > 1. The negative value /3 = —1.3 for cesium implies that it is a low field 
seeker at the wavelength Xl = 157 nm. 

Finally, to see the effects of large laser powers we plot in Figure Hl^a) the quantum 
contrast and in (b) the transmissivity as a function of the light intensity in G3, while 
keeping the other parameters at the values of the marked points of Fig. HJ This case 
corresponds to the solid lines in Fig. |H While the visibility grows with increasing laser 
power, the transmissivity decreases; the same holds for an intensity variation of Gi. In 
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Figure 4. (a) Computed quantum interference contrast and (b) transmissivity when 
tire tlrird grating laser power is varied. The power is expressed in terms of the number 
of absorbed photons. For the solid line in both plots the final signal is given by the 
transmitted neutral particles. The dashed line represents the case of an 'inverted' 
grating G3. Here, the ions are counted directly. The increasing laser power then 
broadens the effective slit width of the third grating pulse. 

this regard, it is notable that the experiment can also be set up such that the count rate 
increases with increasing laser power in the third pulse, by counting the charged clusters 
directly. Even though this 'inverse grating' configuration optimizes the detected cluster 
flux, it again reduces the fringe contrast, as shown by the dashed lines in Figure HI since 
now any increase in laser intensity will broaden the 'effective slit width' in G3. 

3. Theoretical description 

In the following section we recapitulate the description of complex gratings that are both 
ionizing and phase shifting. They are combined to the three-grating interferometer which 
we then model in phase space P ^B51 - HB] to reconstruct the final particle distribution 
and the detected interference signal. Finally, we also incorporate elastic laser light 
scattering. 

3.1. Complex light gratings 

Given a normalized transverse mode profile f{y, z), with J dydz f{y, z) = 1, of the laser 
beam the three-dimensional intensity distribution of the standing light wave is described 
by 

I (x, y, z, t) = 4PL(t)/(y, z) cos' i-j-j , (11) 

where the total absorption is governed by the pulse energy Ei^ = J dtPi^it). There is 
no need to average over the spatial laser beam profile since we assume the clusters to 
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be centered in the laser beam. In the case of a Gaussian laser beam the profile function 
is described by 

f{y,z) = exp . 12 

irWyWz \ '^y "^z J 

The cloud extension is of even lesser importance for pulsed lasers prepared with a flat-top 
spatial intensity profile. 

If the particles are effectively at rest during the nanosecond pulse duration, we can 
assign a complex grating transmission function t^^\x) to each laser pulse k = 1,2,3, 
which transforms a particle's wave function ^(x) — )■ t^''\x)ip{x), according to a 
scattering model |16]. The modulus square < 1 then gives the probability for 

a particle to be transmitted by the grating. The mean number of absorption processes 
per cluster and pulse reads 

= i^±|^t;(o. 0) cos' (H) , „« cos' (H) , (13) 

with the ionization cross section cXabs- The pulse energy E^'^'^ can be varied from pulse to 
pulse. Given the dielectric function e = ei + ie2 of the cluster material the absorption 
cross-section of a spherical sub-wavelength cluster reads [17] 



(Tabs = 47ri?V = ^ / ^0 2^ 2 > 

where R is the cluster radius and q = Sm/AnR^ its mass density. 

The probability for not absorbing any photon during the light pulse is |t^^^(x)p = 
exp(— '^'■'^''(x)). This can be expanded in a Fourier series 

\t^'\^)\'= E 5f(0)exp(^), (15) 

n=—oo ^ ^ 

where the Fourier coefficients are given by Eq. ([6]). 

In addition, the transmission function t^''^ (x) also carries a phase due to the dipole 
energy [2llll6lll8]. It is mediated by the real part of the particle's dipole polarizability 
ttsi(AL) at the laser wavelength, which we use in volume units, a = asi(AL)/47r£:o- For 
sub- wavelength clusters it is given by 

We obtain the phase shift (j)'^^\x) = (f)'^^ cos^(7rx/(i), with 0q'^'' = n^^^ /2/3 and 
This leads to the total grating transmission function 



t'^'\x) = exp ( + 10 W (a;) ) = exp 




2 /7rx\ 
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Its Fourier coefficients are 

= exp (-^ + 1 4 I + I . (19) 




In contrast to the Kapitza-Dirac Talbot-Lau interferometer [TH[T5] where photoabsorp- 
tion led to a contrast-reducing momentum recoil instead of the removal of molecules 
from the beam, both the phase-imprint and the single-photon ionization in G2 are 
independently capable of causing matter-wave interference in our present OTIMA in- 
terferometer. 

3.2. Phase space model 

Starting from an initially bunched ensemble of particles at t = we are now going to 
derive an expression for its final state at t = Ti + T2 using a one-dimensional phase 
space description of its time evolution [T51I211I1S]- The motional density matrix p can 
be rewritten as the Wigner function 

w{x,p) = —— I ds exp (ips/h) {x — ^s\p\x + ^s), (20) 
zvrn, J 

which is a real function of the phase space coordinates x and p. We assume the initial 
Wigner function to be spatially constant on the scale of the grating period d, since 
the incident cluster cloud with spatial extension Aa; is uniformly distributed over many 
periods of the standing wave, Ax ^ d. The momentum dependence is given by the 
one- dimensional marginal distribution 

D{p) = J dpydp,p{p,py,p,), (21) 

where p{p) represents the three-dimensional momentum distribution of the particle 
cloud. The Fourier transform 

D{s) = j dpD{p)exp(-ips/h) , (22) 

which is normalized to -D(O) = J dpD{p) = 1, characterizes the transverse coherence 
of the ensemble. We assume the initial ensemble to be incoherent on the scale of the 
grating period d, i.e. the function D{s) is sharply peaked around s = and nonzero only 
for arguments |s| <^ d. This is equivalent to saying that the momentum distribution 
D{p) is broad and approximately constant on the scale of the grating momentum, i.e. 
Ap ^ h/d. This is the reason why a single grating does not suffice to observe quantum 
interference starting from such an initial state. 

Once subjected to the ffist grating transformation, the initial state wo{x,p) = 
D{p)/Ax undergoes a convolution in phase space Wq{x,p) — j- / dp^T^^^x.p — 
Po)wo{x,po). The grating kernel T^^\x.,p) relates to the transmission function t'^^\x) of 
the fc-th grating via 

TW(x,p) = ^y'd.exp(^^p.)t('^)(a:-0 ■ (23) 
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After Gi the state propagates freely during the time Ti. This corresponds to a shearing 
in phase-space, x — )■ x — pTi/m. The same convolution and propagation transformation 
then applies to the subsequent diffraction due to G2 and the following evolution during 
T2. In order to obtain the final position distribution we integrate the Wigner function 
wti+T2{x,p) over the momentum to get the spatial probability density 

wt^+T2{x) = j dpdpidpoT^^^ (^x - '^,p-pi^ 

xTW (x- P^'^P'^\ p^-po)Dip,). (24) 

Due to the periodicity of the grating kernels f l25]) . wt^j^t-Ss^^ is a d-periodic oscillatory 
function in x. It can be stated in terms of the Talbot-Lau coefficients of the fc-th grating 

00 

(0 = E (^f-n) * [^^^ - 2J)] , (25) 

j=-oo 

which involves the Fourier coefficients from Equation (1191) . We allow for different laser 
powers, i.e. different parameters n^^^ ,(\)^^^ , and following the procedure described in 
we arrive first at an expression for the Talbot-Lau coefficients Bn\0 (Eq- and then 
at the spatial probability density 

1 >^ f2mix\ -fd{nT^+m)\ 



X 



3.3. Resonance approximation 

A clear quantum interference pattern can be observed in the density distribution (126|1 
when the pulses are at least separated on the order of the Talbot time. In this case the 
sharply peaked function D reduces the range of summation indices (n, i) that contribute 
significantly to wti+t2{^)- Only those {n,i) count, which fulfill \nTi + iT2\ ^ T^. For 
Ti,T2 ~ Tt this generally implies a unique integer n for each i [2B\- Here, we restrict 
ourselves to the case, where the delay T2 after the second grating pulse is close to an 
integer multiple of the first pulse separation Ti = T, T2 = NT + r, with |r| <^ T. The 
double summation is then simplified by the resonance approximation 

O (*!Zi±M) . (|I) . (27) 



T / \ T / 



We then arrive at 



- I tf (f ) (I) (^) e.p 

for the position distribution immediately before the third laser pulse. In the resonant 
case, r = 0, we find the ideal Talbot-Lau interference pattern as described by ([2]) for 
N = 1. As noted previously, the Talbot-Lau coefficients of the first grating Bn\o) 
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then reduce to the Fourier components ([6]), and the first pulse serves as a classical mask 
without any phase modulation. This also shows that for an interference pattern to be 
formed, the first grating must not be a pure phase grating. 

This is, however, no longer true for r 7^ 0, where a transient near-field interference 
effect emerges because of the phase modulation at both the first and the second grating. 
This phenomenon was observed in Talbot-Lau interferometry of thermal atoms using 
pure phase gratings [271151]. The maximum time span for such a transient interference 
signal is limited by the width of the initial momentum distribution, |r| < md/Ap. Given 
a typical velocity spread of Av = Ap/m = Im/s, transient interference effects would 
occur for time shifts r of mostly a few hundred nanoseconds with optical phase gratings. 
Absorptive gratings, on the other hand, offer a resonant Talbot-Lau interference effect 
with high contrast, as shown in Section [2751 and discussed below. 

3.4- Fringe shifts and the role of the third grating 



In our previous derivation of fl28l) we assumed the particles to propagate freely between 
the grating pulses. Now, we admit a constant force F = ma along x acting during the 
entire pulse sequence. The propagation in the force field then causes an accelerated 
shearing of the Wigner function 

/ a 2 ^ 

Wt{x,p) = Wq { X 1 — t ,p — mat . (29) 

\ m 2 J 

Our earlier calculation can be redone for all pulses and reveals a shift x ^ x + 6x of the 

density pattern fl2S]) at time (A^ + 1)T, given by 

Sx = ~N{N +1)T^. (30) 

Since all three grating laser beams are retro-refiected by the same fixed mirror, we may 
now exploit the fringe shift, Eq. fl5U]) . to scan the interference pattern by tuning the time 
separation. This will shift the fringe coordinate xs and modify the Fourier components 
of the detection signal after the third pulse. 



|2 



S{xs) = J dxw(^N+i)T{x + xs)\t^ '{x)\ = J dxw(N+i)T{x)\t^ '{x - Xs)\ . (31) 

In practice, the Earth's gravitational acceleration a = g provides a highly homogeneous 
and constant force that can be easily used to shift the fringe pattern if the x-direction 
is chosen to have a vertical component. Alternatively, one may also use electrostatic 
fields. 

When the signal is recorded as a function of T it will no longer be strictly periodic. 
In practice, however, the time variation d/N{N + l)gT required to shift the pattern 
by about one grating period d can be made small compared to the Talbot time Tt by 
increasing either the pulse separation time T or the mass of the particles. The sinusoidal 
visibility then remains a good measure. After inserting fl28|) into fl3T]) . we obtain the 
final signal 
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Its fringe contrast is 

B% (r/Tr) {{NT + t)/T^) B^}{0) ^ f dr 



(33) 



5«(0)5f (0)i?J^)(0) \Tt, 

This shows that it is favorable to keep the pulse delays Ti and T2 equal, i.e. to set = 1: 
Any larger N increases the index of the Talbot-Lau coefficients ([3]) and thus reduces 
the maximum fringe contrast. In addition, the simulations of Section 12.51 show that a 
high visibility is realized with ionizing gratings in the resonant case r = 0. Equation 
032]) suggests that this deteriorates rapidly with increasing |r| because of the sharply 
peaked Fourier transform D{s) of the momentum distribution. We therefore consider 
the resonant case r = with equal pulse timing, = 1, in the following. The signal 
expression (15^ and the visibility (155]) then reduce to Equations ([7]) and respectively. 

In an alternative 'inverse configuration' of G3, where the ions are counted instead of 
the neutral clusters, we have to replace |t''^^(a;)p by l — \t^^\x)\'^ . The grating coefficients 
in ^ then change from Eq. (E} into Bi?\o) = 6^,0 - exp{-n^Q^ /2)In (^-n[,^V2)- 

4. Challenges and limitations 
4-1- Rayleigh scattering in the grating 

So far the influence of the optical gratings was treated as a coherent transformation that 
modulates the cluster wave function through the complex transmission function ( IT8|) . 
However, in addition to the absorption of photons and their virtual scattering within the 
laser wave the diffracted nanoparticles may also give rise to Rayleigh scattering, where 
the photon is re-emitted elastically. This process causes momentum diffusion P|H^ and 
thus deteriorates the interference pattern. 

For short laser pulse durations we can treat the Rayleigh scattering independently of 
the coherent grating transformations. We model it by a Lindblad-type master equation 
for the density matrix [501 - I52] 

dtp = 7r f du N{u) cos ( ^lx) e-'^'L«-Xpe''=L"^'^ cos ( ^lx) - ^ { p, cos^ ( A;lx) } 
iJ\u\=i 2 

The position operator is denoted by x. The equation describes the elastic scattering 
of a single photon from the diffracting laser fleld into a random direction u. The 
corresponding momentum recoil is represented by the operators for the standing- 
wave mode cos/clx and the plane wave modes exp (i^LMajx). The number of recoils 
is determined by the total scattering rate 7r at the antinodes of the light fleld, and 
N{u) gives the distribution of scattering directions, as described by the dipole pattern 
A^('u) = 3 sin^ -i^/STT. Here 'd is the angle between u and the polarization vector of the 
laser fleld. 

We can solve the master equation flM|l in the position representation. The density 
matrix then transforms according to {x\p\x') — )■ ri{x,x'){x\p\x') , with the decoherence 



(34) 
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function 

ri{x, x') = exp 



^ ( 2 / du A^(u)e*''^""^^' cos ki,x cos kj^x' - cos^ ki,x - cos^ ki,x' 
, 2 \ J\u\=i 



.(35) 



It is trace-preserving, ri{x,x) = 1, and it only acts on the spatial coherences. The mean 
number of scattered photons at the antinodes of the kth grating is related to the number 
of absorbed photons via the ratio of the associated cross sections, n^^ = n'^Q ^'an/aahs- In 
phase space representation, decoherence due to elastic scattering is then represented 
by the integral kernel 



2nh 

oo 

E 

n=— oo 



ds exp 
/ 27rinx\ 



1 , (k) f S S 

-ps + n\^'r] l^x- -,x + - 



exp j y d. exp (^ip.j Rl^^ (^) 



(36) 



which has to be applied at each grating pulse k in addition to the coherent grating 
kernel (l23l) . Again, the periodicity of the kernel in x allows us to expand it as a Fourier 
series, where the closed expression for the Fourier coefficients 



exp 



X /„ 



3n 



{k) 

R 



3n 



(k) 



cos TT^ ( jo (ttO - 

Ji (vrO 



Jo (vrO - 



2 , 
- cos vrt 

3 ^ 



(37) 



follows from a lengthy calculation involving combined integrals over trigonometric 
expressions and Bessel functions. The symbol j denotes the spherical Bessel function of 
the first kind. Combining the coherent kernel fl23|) and the scattering kernel fl36|) results 
in the modified Talbot-Lau coefficients 



(38) 



for each grating. Since Bn\o) = Bn\o) the Rayleigh scattering acts as a purely 
decohering agent and does not infiuence the periodic masking by the gratings. This 
means, in particular, that only the treatment of the second grating needs to be modified. 

With growing cluster size the role of elastic light scattering becomes increasingly 
important. Since the scattered photons have the same short wavelength as the grating 
light, they can decohere the previously delocalized matter waves. We may disregard 
this effect only as long as the number of scattering events is small compared to the 
number of absorption processes n^^K Their ratio is given by their total cross sections 



n 



(2) 
R 



n, 



(2) 



Cabs 



2 (51 - 1)^ + el 



9 



£2 



{kLRf 



2 m 



£2 



(39) 



and it is proportional to the mass, assuming that the dielectric function e = Si + 182 
and the density q are constant. The influence of Rayleigh scattering on the interference 
contrast is shown in Fig. El where the expected sinusoidal fringe visibility is plotted as a 
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pulse delay time T/T^ 



Figure 5. Calculated fringe visibility including Rayleigh scattering of laser photons as 

(2) 

function of T/Tt for a fixed Uq ' and gold clusters of different masses. The solid and the 
dashed lines correspond to m = 10* amu and m = 10^ amu, respectively. The dotted 
line reproduces the quantum contrast curve of Fig. [^Ja) without Rayleigh scattering, 
which is a good approximation for cluster masses < 10* amu. 

function of tlie pulse delay around the Talbot time, both for gold clusters of m = 10^ amu 
(solid line) and for m = 10^ amu (dashed line). The power of the UV laser gratings is 
adjusted such that the absorption parameters are kept constant at rS^^ = n^^^ = n^^^ = 8 
(as in Section [275]) . While the solid line still resembles the predicted visibility without 
Rayleigh scattering (dotted line), the overall contrast of the dashed curve is somewhat 
suppressed due to the increased scattering cross-section of = 0.9(Tabs at m = 10^ amu. 
This corresponds to a gold cluster of 55 nm diameter, a size where the point particle 
approximation ceases to be valid. It also exceeds by two to three orders of magnitude 
the mass limit where the experiments have to be modified in order to compensate for 
the influence of gravity (see Sec. 14. 3p . 

4.2. Collisional decoherence and thermal emission 

Another concern in interferometry with massive particles is their uncontrolled 
interaction with the environmental degrees of freedom. The most prominent causes 
of decoherence are collisions with residual gas particles and the emission or scattering 
of thermal radiation |53l|5l]. We find that the unperturbed interference of a 10^ amu 
gold cluster after a delocalization time of 2Tt ~ 30 ms requires that the background gas 
pressure is kept at p < 10~^ mbar and the internal cluster temperature at T < 1000 K. 
Both values can be readily achieved in an experiment. For this estimate we chose N2 as 
the background gas and use the London dispersion formula to compute its interaction 
strength with the gold clusters, assuming the ionization potential of N2 to be 15.6 eV and 
its static polarizability oatj = 47r£:o x 1.74A^. The dielectric function of gold at thermal 
photon wavelengths is here assumed to follow Drude's model with a work function of 
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W = 5.4 eV, a plasmon frequency of up = 1.3 x 10^^ Hz and a resonance width of 
Fp = 1.1 X 10^^ Hz. The static polarizabihty of gold is taken to be that of an ideal metal 
cluster. 

4- 3. Fringe averaging on Earth 

As seen in Sec. 13. 4|. gravity can cause a fringe shift if the direction of the force is parallel 
to the diffraction grating vector. For usual interferometers this phase shift is dispersive, 
i.e. larger for clusters of lower velocity and longer falling time. Even a small width of 
the cluster velocity distribution can then cause a rapid fringe averaging in the plane 
of G3. This problem can be completely eliminated in an OTIMA interferometer where 
all clusters are exposed to the gravitational field for exactly the same amount of time. 
All phase shifts are therefore equal and one can still expect high contrast interference, 
independent of the grating orientation in space. 

What remains is the issue of the classical falling distance. At 10^ amu a cluster 
needs a passage time of about 30 ms, assuming F2 laser gratings with a period of 78.5 nm. 
During this time it will fall by about 4.6 mm. At m = 10'' amu and a time of 300 ms 
the free fall distance reaches already half a meter. This illustrates that the force of 
gravity will have to be compensated at some point by electrical or magnetic fields or to 
be eliminated in a microgravity environment. 

In contrast to gravity, the Coriolis force Fq = 2mv x Q depends on the cluster 
velocity both in modulus and direction. It cannot be compensated by clever timing. 
However, OTIMA offers still an advantage over a fixed-length interferometer because it 
allows us to align the mirror orientation normal to Q. The resulting Coriolis force is 
then intrinsically orthogonal to the grating vector d, such that the acceleration does not 
contribute to the fringe shift. 

With growing masses and Talbot times the interferometer becomes more sensitive 
to vibrations. For low-frequency vibrations, the total fringe shift depends on the position 
change Ad of each individual grating according to Ad = Ad^^^ — 2Ad^'^^ + Ad^^\ The 
single mirror design of the proposed OTIMA interferometer (Figure [T]) ensures that 
many conceivable contributions drop out, such as quasi-static mirror shifts and mirror 
tilts. All dynamical mechanical effects have to be suppressed on the level of about 

5- 10 nm, which can be done using established spring suspension systems. 

5. Experimental considerations 

Having described the general idea, we now discuss the compatibility of existing beam 
sources, diffraction elements and detectors with our suggested scheme and we start with 
the availability of single-photon ionization (SPI) gratings [2Z] for various materials. 
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5.1. Grating requirements 

For our present proposal we consider only particles, for which efficient single-photon 
ionization (SPI) has already been observed [55I457] . such as metal clusters. The cluster 
ionization energy usually exceeds the work function of the bulk by less than 1-2 eV. It 
decreases with the particle radius like i^ion — W + 0.426^ /AtteoR and approaches the 
bulk value at 20 to 100 atoms per cluster [52]. Dozens of materials are therefore suitable 
candidates for SPI diffraction and detection, provided that the photon energy is higher 
than 5eV, i.e. A < 250 nm. 

Although SPI yields are not quantitatively known for all materials, the process has 
been shown to be dominant in several systems, also in competition with fragmentation 
and photoemission. This applies to semiconductor nanoparticles [58] as well as to metal 
clusters [59]. In our simulations (Figs. El |3l S]) we have assumed a yield of ?7spi — 1. 

In the following we distinguish interferometry with particles below 20,000 amu and 
from 10^ to 10^ amu, since they may require different implementations of the same idea. 

We envisage a de Broglie wavelength of XdB — 5 x 10~^^m, to fit the entire 
interferometer on the same two-inch mirror. This is compatible with a velocity of 40m/s 
for clusters with m < 20, 000 amu. This layout is driven by the desire to have a short 
and ultra-stable base for all laser gratings, which is important for high signals and with 
regard to a high common-mode rejection for phase fluctuations related to vibrations and 
rotations of the mirror mount. 

For this setting, a F2 excimer laser is well adapted. It emits a few hundred VUV 
light pulses per second with an energy of more than 1 mJ and a duration of less than 
5 ns. Retroreflection of the laser light generates a diffraction grating with a period of 
d = 78.5 nm. The F2 laser wavelength is fixed at A = 157.63094(10) nm with a line 
width of less than 1 pm [HO]- The gain proffie of fluorine allows, in principle, a five times 
weaker line at A = 157.52433(10) nm, as well. The longitudinal coherence length of more 
than 1 cm, guarantees a sufficient grating periodicity up to a distance of 2 mm from the 
retro-refiecting mirror. This covers the entire width of the cluster beam and about 10^ 
grating periods. 

The transverse coherence of excimer lasers is limited and often determined by 
diffraction of the light beam at the laser outcoupler window. Following the van Cittert 
Zernike theorem [61] and also recent measurements the transverse coherence is estimated 
to be 80 /im in a meter behind the laser window [62] . The homogeneity of the standing 
light wave is dominated by the quality of the mirror surface. A flatness of better than 
d/10 is available for VUV wavelengths. Given this, even the limited transverse coherence 
is still sufficient for the formation of a homogeneous diffraction grating, provided that 
the laser beam divergence is limited to about 1 mrad. 

All this suggests that a F2 laser is the ideal basis for building a compact and 
rugged OTIMA interferometer which allows one to address a large range of materials. 
This flexibility is bought, however, at the price of working with vacuum ultraviolet laser 
light, which requires purged beam lines and specialized optics. 
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The short laser pulse width is an advantage for precision measurements but also 
entails some further consequences related to the fact that the standing light wave 
needs some time to form and that it also experiences an overall temporal amplitude 
modulation. The finite speed of light is the reason why the clusters will interact with a 
running wave during a fraction of the laser pulse duration. For a cluster beam confined 
to within 1mm from the mirror surface, this running wave is limited to 0.5% of the 
total pulse energy. Also the low reflectivity of even the best mirrors which are available 
for 157 nm {R > 98 %) causes the cluster beam to be exposed to a running light wave of 
2 % the total pulse energy. This slightly affects the final count rate but not the visibility, 
since only the neutral particles reach the detector. 

Amplitude variations of the laser light field, realistically on the order of 3 to 
5%, modulate the population of different diffraction orders but they cannot destroy 
the interference pattern, since the fringe positions are solely determined by the laser 
wavelength. Also the pointing instability of the grating lasers is negligible if one can 
ensure that all detected particles have interacted with all three gratings. A critical point 
is the accurate timing, i.e a low jitter, of the laser gratings. With existing technologies 
a timing accuracy of 1 ns can be achieved, which is sufficient for the purpose. 

Although the Talbot time = (Pm/h is intrinsically independent of the particle's 
velocity, clusters of different speed will travel different distances between the gratings of 
the same interferometer. Therefore, the cluster beam must be sliced into packets (Figure 
d]) such that all particles of a batch always interact with the same spatially extended laser 
pulse, independent of their velocity. The combination of a chopped source with a pulsed 
TOF-MS detector will, in practice, select a velocity band of Au/w ~ 2% in the proposed 
experiments up to 10^ amu. A selection of de Broglie wavelengths is then achieved 
in combination with the detecting mass spectrometer. In principle, an interferogram 
may then even be recorded without scanning either velocity, time or grating position, 
provided the source emits a sufficiently dense set of cluster masses. 

In the range of 10^ amu the particle cloud will already expand quite significantly 
during a total interferometer time of 30 to 60 ms, even if we assume a cluster velocity of 
Im/s, which is a challenge in itself. We thus require a laser with a spectral coherence 
of Ai/ < 0.1 cm~^, as for instance provided by a grating stabilized ArF laser at 193 nm, 
the fifth harmonic of a seeded Nd:YAG laser at 213 nm, a frequency doubled dye lasers 
or an optical parametric oscillator in the wavelength range between 210 and 250 nm. 
Since the absorption cross sections of metal clusters grow linearly with the mass for 
a given atomic species, their higher values will then compensate for the intensity loss 
when the UV laser beam must be expanded to cover the growing cluster cloud. Finally, 
in the high-mass range the temporal envelope of the laser beams grows in relevance. 
The expanded cloud will now see a running wave for up to several percent of the total 
pulse width. This time may be reduced by making a standing light wave from counter 
propagating laser beams. 
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5.2. Source properties 

On the source side, different options are at hand for different mass ranges: Interestingly, 
even an effusive beam source may match the requirements of an OTIMA interferometer. 
The most probable velocity for m = 10^ amu at T=500K is v^p ~ 3m/s and recent 
experiments have shown that there is still hope for the thermal volatilization of highly 
massive but chemically tailored organic molecules [53] . 

However, for the purpose of our present proposal we restrict ourselves to a cluster 
condensation source, which first volatilizes the atoms and then recondenses them in a 
cold noble gas stream. The initial atomization may be done thermally, for some materials 
such as alkali atoms [55] and fuUerenes [61], or using a magnetron sputter technique for 
a large range of metal clusters with higher melting points [65j . A gas aggregation source 
ejects a broad distribution of clusters, ranging from a few up to several ten thousand 
atoms per cluster with a brilliance of up to 10 mg srad/s. 

Own preliminary studies show that under conditions comparable to OTIMA we 
can count about 100 cluster ions on the TOF-MS MCP in a time window of 200 fis. To 
achieve an adequate visibility each laser grating has to reduce the cluster beam by a 
factor of three. After the interferometer we thus retain 4% of the initial beam, i.e. 600 
ions per second with virtually no background during the measurement period, which is 
sufficient for the purpose of interferometry. 

In the setup of Figure [1] the cluster velocity is determined by the flow of the cold 
carrier gas, which is usually a mixture of Helium and Argon. The most probable cluster 
velocity in a nitrogen cooled setup amounts to about v^p = 300 m/s. This can be 
lowered by a cryogenic cooling stage, using a neon buffer gas at 30 K [66J to achieve 
a cluster velocity in the range of less than 100 m/s, if we include the relative velocity 
slip of massive clusters in light noble gases ^7\. From this distribution one may still 
select the low- velocity tail to get to particles in the range of 40 to 50 m/s. Moreover, 
experiments aimed at surface cleaning also showed that laser acoustic desorption is 
capable of generating nanoparticle beams with speeds in the range of 10 to 50 m/s |68j . 

Many complementary methods are still conceivable. Of particular interest are 
those techniques that focus on the improvement of cooling for freely suspended 
mesoscopic particles [69H72]. Although most published proposals are still best adapted 
to small source volumes and a few particles, a rapidly growing community is currently 
investigating methods for preparing mesoscopic isolated particles, which may eventually 
also be coupled to our proposed interferometer scheme. All this suggests that there are 
ample possibilities for further progress to masses up to 10^ amu. 

6. Conclusions and Outlook 

The proposed OTIMA-interferometer is universal in the sense that a single experimental 
setup will enable experiments with a large class of nanoparticles, ranging from 
single atoms to metal clusters, semiconductor nanocrystals, and possibly even some 
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biomolecular complexes, as long as ionization or neutralization is the dominant response 
to the absorption of a single photon. This covers a mass range from 6 and 7amu for 
lithium isotopes to 10^ amu for cold giant clusters. 

The setup is promising for testing the foundations of wave mechanics, and its 
intrinsic features, such as high common mode rejection, non-dispersiveness and good 
control over the grating pulse times, are well-adapted to interference-assisted cluster 
metrology. The OTIMA concept could thus enable new insights into cluster physics, 
including coUisional and absolute photoabsorption cross sections, electric and magnetic 
properties etc. 

Our calculations show that a high interference contrast may be expected for a wide 
range of masses under otherwise very similar conditions. This could become relevant in 
explorations of the equivalence principle as well as in experiments that probe currently 
speculative fundamental limitations of the Schrodinger equation [TSHZH]. 

In contrast to many other experiments where source purity is a value in itself, the 
OTIMA interferometer can profit from the broad distribution of particles ejected by 
a magnetron clusters source [HS]. Performing experiments with different cluster types 
under otherwise identical conditions provides a high degree of 'common mode noise 
rejection' and a good way of identifying systematic phase shifts under variations of 
mass, baryon composition, geometrical structure or other cluster properties. 
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